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Abstract 

> 

CO ' We consider the Casimir energy of a thick dielectric-diamagnetic 

shell under a uniform velocity light condition, as a function of the radii 
qq , and the permeabilities. We show that there is a range of parameters 

in which the stress on the outer shell is inward, and a range where 
the stress on the outer shell is outward. We examine the possibility of 
obtaining an energetically stable configuration of a thick shell made 
of a material with a fixed volume. 



1 Introduction 

It is well known that the fluctuations of electromagnetic fields in vacuum or 
in material media depend on the boundary conditions imposed on the fields. 
This dependence gives rise to forces which are known as Casimir forces, acting 
on the boundaries. The best known example for such forces is the attractive 
force experienced by parallel conducting plates in vacuum . Casimir forces 
between similar, disjoint objects such as two conducting or dielectric bodies 
are known in most cases to be attractive [0 and are sometimes viewed as the 
macroscopic consequence of Van der Waals and Casimir-Polder attraction 
between molecules. 
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However, Boyer || showed that the zero point electromagnetic pressure 
on a conducting shell is directed outward^ In this paper we address the 
question: can there exist a compact ball for which the Casimir forces would 
not expand the ball to infinity? We look for such behavior in a simple model. 

In view of the dominance of the Casimir forces at the nanometer scale, 
where the attractive force could lead to restrictive limits on nanodevices, 
the study of repelling Casimir forces is of increasing interest. Indeed, Boyer, 
following Casimir's suggestion, studied inter plane Casimir force with one 
plate a perfect conductor while the other is infinitely permeable. He showed 
that, in this case the plates repel ||. This problem was reconsidered since 
in §B. 

In addition to the Casimir effect for a conducting spherical shell, various 
cases of material balls where considered in the literature. The case of a ball 
made of a dielectric material was considered by many authors and exhibits 



strong dependence on cutoff parameters [0, |9|, [I(J [IT], 0. The case of 
a dielectric-diamagnetic ball has also been extensively studied, especially 
under the condition of e/i = 1, which will be referred to as the uvl (uniform 
velocity of light) condition, since its introduction by Brevik and Kolbenstvedt 
[p~3| , |I3| , [I5| , |l6f| . A medium with the uvl property has in many cases cutoff- 
independent values for the Casimir energy (see also 18, [H| for the Casimir 
energy of a cylinder with uvl). A heuristic argument for this statement 
goes as follows: the zero point energy of the electromagnetic field in a uvl 
medium is a sum over the eigenfrequencies of the system, and behaves for 
high frequencies as Yli^i ~ Yli c ki, where the factor c is common to all the 
media and the geometric information on the system enters via the allowed 
k's. This expression in the limit of high k's behaves exactly as the vacuum 
energy with e/i = c~ 2 and the sum becomes regularizable by subtraction 
of the vacuum energy. In cases where this condition is not fulfilled there 
are problems of UV divergent terms which are proportional to the volume 
in which the velocity of light is different from the velocity of light in the 
background. 

We will use the uvl condition for our thick shell. 

In all of the above mentioned cases, namely, the conducting sphere, di- 
electric ball or a dielectric-diamagnetic ball with uvl, the resultant pressure 
was found to be repulsive. 

In order to obtain diverse behaviors we mix the case of disjoint bodies 



2 Assuming that zero point forces are in general attractive, Casimir considered a semi- 
classical toy-model for the electron in which the coulomb self-repulsion is balanced by a 
Casimir type attraction. One of the consequences of Boyer's calculation is that since the 
pressure on a conducting sphere is outward it cannot balance the Coulomb repulsion in 
Casimir's model. 
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Figure 1: A thick dielectric-diamagnetic shell. 



(where there is usually attraction) and the dielectric ball scenario as follows: 
We consider the case of a thick shell (fig. [I]), with three permeability regimes 
(inner, middle and outer) fij, fin, \iui, under a uvl condition. In this case 
there are two competing effects, i.e.: interaction between the inner and outer 
boundaries, and the repulsive pressure experienced by each boundary. 

Using a formula derived in |2l| we show in sections ^ and |3] that for a 



dilute medium the energy of the thick shell is 

rn / m 5 2 r 3 (r 2 -5R 2 ) 2 
E c (r, R) = — K R + K r K R — — — — + K r - 



USttR R " r H A{r - Rf{r + Rfn " r 128irr 

Where r and R are the inner and outer radii respectively and the parameters 
n r and are defined in section |3|. From this expression one can easily obtain 
the energy of a single ball by taking the limit of R to infinity, and also the 
energy of two parallel infinite dielectrics by taking both of the radii to infinity 
while keeping a finite distance d = R — r. 

Next, in sections |] and [5] we investigate this expression for the energy 
and show that there is a range of parameters such that the force on the outer 
shell is attractive. 



2 Energy of the electromagnetic fluctuations 



In this section we briefly review the Green's function method for calculation 
of Casimir energies (see [^TJ and [^] for details). To calculate the Casimir 
energy of a medium under the uvl condition, we use the perturbative tech- 
nique suggested in p0[. There, the Born series for the correlation function 



Dik(u>; r, r') of the electromagnetic fields in material medium is presented. 
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The correlation function Dik is defined by 

POO 

A fc (^;r,r')= / e^Dg(t;r,r')dt (1) 
Jo 



where 



nR/ . , . f!<ii(r,ti)4(r'i2)-4(r't2)A(r,ii)> h<t 2 
V ik {t 2 -ti,r,r ) - | Q otherwise 



is the retarded Green's function. Throughout we use the gauge A = 0, so 
that the indices i, k range 1, 2, 3 and D is a 3 x 3 matrix. 

The correlation function D is known |H] to be the Green's function for 
the equation (in units where H — c—1) 

([V]/i -1 [V] - u 2 e)D = -I (2) 

where [a] stands for a matrix whose elements are [a]ik = tijk a j and I is 
the identity operator, which, in coordinate space is just the 3x3 identity 
matrix times a delta function. In eq. e is the permittivity and \x the 
permeability of the medium. In the following we assume that /z(r) and e(r) 
are scalar functions (of course, in the general case both /i and e are tensors). 
This equation can be also written in the form: 

([V] 2 - [Vlog/x][V] - u 2 efj,)D = -//I (3) 

The correlation function of electromagnetic fields in vacuum, where e = 
/i = 1, will be denoted D . It is the inverse of the operator ([V] 2 — uj 2 ) and 
is well known. It is given by the formula: 

D (u; r, r') = - (l + \ V ® v) - , 1 - exp^l r ~ r 'D, (4) 
\ uj z J 4-7r|r — r | 

We now wish to use the known Dq to express D via a Born series. Define: 

V = u 2 I(l - fie) (5) 

and 



Q = -[Vlog(/i)][V] (6) 
It follows from (|3|) that as an operator D satisfies: 

D = ((I - (V + Q)D )D - 1 )- 1 fi = D (I -(V+ Q)-D ) - V> (7) 
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Thus D is given by the following formal series: 



D = D of i + D {V + Q)D of jt + D (V + Q)D {V + Q)D fi + .... (8) 

We now impose the "uniform velocity of light" condition in the medium 
by setting epi = I. This eliminates the V terms and we are left with the 
expansion: 

D = D /i + D QD q ij, + D QD QD ii + .... (9) 

The correlation function D can be used to calculate various properties of 
the electromagnetic field. We use it to calculate the energy density of the 
field by using the relations: 

(EiE k ) u = u^AAk)^ 

and 

(BiB k ) u = V x (AiA h ) u x V, 

Where the round brackets stand for the Fourier transform with respect to 
time of the correlation function | < Aj(r, t)Ak(r / , 0) + A k {v' , 0)Aj(r, t) > . 

At a temperature T = -s this correlation function of the fields At is related 
to the retarded Green's function D by the fluctuation dissipation theorem 



(A(r)^(r')). = -z47rcoth(^)(A,(r,r') - ^(r',r)) (10) 
Thus, at zero temperature, the energy density of the electromagnetic field is: 

^lim r /^ r ImTr^ 2 e(r)D(r,r / ) + x D x Vjda;, 



where we have chosen to neglect the dependence of the permeability \i and 
permittivity e on the frequency u. Inserting the expansion @ in (|TT]) one can 
obtain a series for the energy density of fluctuations of the electromagnetic 



field. It was shown in [2(J that the first contribution to the Casimir energy 
density (i.e., energy of the electromagnetic field in the presence of exter- 
nal conditions minus the energy density of the electromagnetic field without 
them) comes from the term D QD QD fi in (0) Q For a dilute medium this 
term gives the dominant contribution to the Casimir energy. 

3 To see this, note that the Do term is just the contribution of a homogeneous vacuum 
which we subtract. The terms of the form DoQDoQfj, can be shown to cancel between the 
electric and magnetic correlation functions. 
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3 Casimir energy of a thick shell 

The Casimir energy per da; is obtained via eq. (|ll|) by the integration of 
p(r,u) — (contribution of D ) over space. This yields the following general 
formula for the density of Casimir energy of a medium with a radially sym- 
metric permeability \i |2(| assuming the uvl property, and diluteness: 



Where: 



p T ( 2 )( w ) = _Wcoth(^)Im [°° ds I (s) (log fi(s))' (12) 

2 Jo 



1(a) = ^dudv { - ^ S)) g 2 (u)(l0g»(v))' (13) 

and the integration domain T is such that u, v and s can form a triangle. 
The function g 2 is given by: 

We consider a shell of thickness R — r with inner radius r, and an outer 
radius R. The permeability of the shell is \xu (Fig. |l|), it is imbedded in a 
medium with permeability \xm and its core has permeability fij. We define 
K r = log — and % = log iiLL . 

1 ° MIX Vin 

We now calculate the Casimir energy of the thick shell using eq. QI^D- 
In our case \i is a sum of two radial step functions, and fjf(s) is just a pair 
of delta functions at s = r and s = R. Thus the integration over s in (jl2|) 
becomes immediate and yields the energy density 

p T ( 2 ) ( u ) = -Attu 2 coth(^)Im (n r I(r) + krI(R)) (15) 



The Ps in can be explicitly calculated: 



; 2r , . (m 4 - (R 2 - r 2 ) 2 ) 

./(/•) = k,./ y 92 ^ + Kr J r ^ ^ ( " ) 



2m 



and: 



So that 



I(R) = k r / — g 2 (u) + K r / ^ K - ! -^g 2 (u) (17) 

J0 1 JR-r Zu 



= -4 7 ro; 2 coth(^)Im(/ t 2 / 2r ^ 2 (n) + (18) 

2 r2R u 3 / \ , rii+r (« 4 -(R 2 -r 2 ) 2 ) / n\ 
fl Jo T&W + J R _ r 1 2» 92(U)) 



K 
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We can identify the first two terms as the densities of Casimir energy per da; 
for balls of radii r and R. The total Casimir energy can now be obtained 
by integration of the density flTgp over the frequencies u using ( GE| ) and (E7p. 
The result is: 

„ , „ N 5 „ r 3 (r 2 -5R 2 ) , 5 

GV ' ; 128ti\R K w 4(r - i?) 3 (r + i?)% r 128^r v ; 

This is our final expression for the energy of a thick shell, where the assump- 
tion of diluteness implies |« r |, |«r| << 1. 

Let us check that this result coincides with the known result for a single 
sphere when R goes to infinity. To do so we expand K r in terms of the 
diluteness parameter £ = which is commonly used in the literature: 

K r = log = log = -2£ - 2^ + ... (20) 
fin l + £ o 



Substituting this expansion in fllTf ) and taking the limit R — ► oo we imme- 
diately regain the usual result for the energy of a dielectric-diamagnetic ball 
|16| ) p2]| , namely: E c (r, oo) = + 0(£ 4 ). This energy yields an outward 
pressure on the ball. 

Note also that the Casimir energy of two parallel dielectrics per unit area 
can be obtained from (|P9"|), by taking the limit of large radii. To see this we 
keep d = R — r finite while taking the limit r, R — > oo, and divide by the 
surface thus obtaining;: 

*«> - (21) 

To illustrate the broad range of behaviors that are possible from an ex- 
pression for the energy such as (0), we study in more detail two cases: 

1. K r — — Kr. This happens when the inner and outer material are the same. 
E.g. one can imagine a thick material shell in vacuum. 

2. K r = kr. In this case the ration of magnetic constants between the inner 
and middle materials is the same as the ratio between the middle and outer 
material. 



4 Stress 

Now we use the expression for the energy ( |T9"D in order to evaluate the Casimir 
stress on the shells. The pressure on the outer shell is given by 

Id 1 /r 3 R(r 2 + 5R 2 )K r K R 5k 2 r \ 

R ~ ~4^dR C ~4^R^\ (r-R)\r + RY + V28R 2 ) ( ' 
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Figure 2: Pressure on outer shell as function of the radii ratio. The Casimir 
force changes direction at — ~ 0.19 



While the pressure on the inner shell is given by 

l _^_ F 1 / r 2 (r 4 -10r 2 i? 2 -15i? 4 )^^ 5k 2 \ 

An 2 r 2 dr C An 2 r 2 \ 4(r — i?) 4 (r + R) 4 128r 2 J [ ' 

One may recognize the self Casimir force acting on the inner and outer shells 
in the terms independent of r or R respectively. We now turn to examine 
two cases: 



4.1 K r = — K.R 

The pressure on the outer shell is 

.2 



T 



/ 5 _ r"R(r 2 + 5R 2 ) x 
! R 2 Vl28f? 2 (r- RY(r + RYJ 1 ' 



Att 2 R 2 \ 128R 2 (r - R)\r + R) 

The sign of this expression depends only on the ratio c = The pressure 
can be written in the form: 

_ k 2 „ / 5 c 3 (c 2 + 5) \ . . 

1^(128 - (C -W+ o<c<i < 25 > 

For a constant outer radius R, the behavior of the pressure as a function 
of the ratio c is indicated in Fig.0. One can see that at a radii ratio of ap- 
proximately 0.19 the Casimir pressure suddenly changes sign and becomes an 
inward pressure. However, the pressure on the inner shell is always directed 
outward so that the inner and outer regions attract. This attraction may be 
balanced by adding the compression resistance of the middle medium, or by 
adding the volume dependence of /i. In the case that r < 0.19i?, the pressure 
on the inner shell is larger than the pressure on the outer shell, so that one 
can imagine the two shells growing until they reach the 0.19 ratio and then 
the outer shell will start contracting. 
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This result is not too surprising. If one considers two conducting shells 
which are very close, then there will be attraction between the shells, of the 
order of magnitude of attraction between conducting plates. This attraction 
will lose its dominance once the radii become far in magnitude and then each 
shell will experience its own outward Casimir pressure. However, in order 
to know if the outer radius stays finite or goes to infinity, it is necessary to 
introduce the law by which the smaller radius changes as we change the outer 
radius. An example for such a situation is given in the next section. [| 

4.2 K r = kr 

In this case the pressure on the outer shell is outward for all R. The inner 
shell will be subject to an outward pressure as long as R > 3.46r. For 
R < 3.46r the pressure will be an inward one on the inner shell. We return 
to this situation in the next section. 

5 A thick shell with a fixed volume 

In the previous section we calculated the pressure assuming that the volume 
of the inner ball is fixed. Another interesting model is to consider the volume 
of the thick shell itself as constant. Thus we can assume v = ^-(R 3 — r 3 ) 
remains constant, and look for a minimum of the energy under this condition. 

In this the outer shell expands, the distance between the inner 

and outer boundaries decreases. So that if there is attraction between the 
shells, it will be energetically favorable for the shell to expand to infinity 
gaining from the energy of interaction between the boundaries as well as 
from the tendency of each separate boundary to grow. This is indeed what 
happens in the case where K r — — Kr. 

However, for the second possibility, namely K r = kr there is repulsion 
between the boundaries. This protects the shell from growing to infinity 
and a stable minimum of the potential can be found. A typical example is 
illustrated in fig. |3|, where the minimum of the potential is for R ~ 1.01, 
r ~ 0.31 assuming the volume of the substance in the middle region is kept 
constant at 

Another possibility is to keep the distance between the shells constant 
(imagine the inner and outer boundaries are attached by means of stiff rods 

4 In this section we considered n r = — kr but the qualitative results will hold whenever 
K r and kr have different signs, as long as they are both small and of the same order of 
magnitude, so as to satisfy the validity of eq. (p"2|). 
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. 18 




Figure 3: The potential of a thick shell with a fixed volume as function of R, 
for n r = kr. There is a clear minimum at R = 1.01. 



of constant length). Qualitatively the results are then similar to the results 
discussed above for the cases K r = — kr and K r = kr. 



6 Discussion 

One of the most intriguing aspects of the Casimir force is its sign: this force 
was shown to exhibit different behavior in different problems. It seems that 
the sign of the force depends on the balancing of several effects, such as a 
tendency to minimize the average curvature of the boundaries, and on the 
other hand the interaction between different patches of the boundary. In the 
calculations here, which are confined to the case of uniform velocity of light, 
we have shown how these effects can work with or against each other in an 
explicit way resulting in a new range of behaviors. We obtained a general 
expression for the energy as function of the radii and permeabilities in the 
limit of dilute media. In particular, for fixed distance between the shells, in 
the limit of radii approaching to infinity, we regain, in essence, the standard 
expression for the Casimir energy between parallel dielectric media. However, 
now the sign of the interaction energy in (|T9|) is determined by the relative 
size of their permeabilities: for a shell enclosed between two vacuua we get 
attraction between the boundaries, as it should be for the standard Casimir 
two plates case. Repulsion is obtained when the permeability of the shell itself 
is between the permeability of the inner substance and the permeability of 
the outer substance (i.e. Hi < fin < nm or /i/ > fin > /im)Q. 

5 The repulsion of plates in such a case can be confirmed by a direct calculation in the 
usual geometry of two infinite dielectric-diamagnetic media seperated by vacuum. 
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7 Appendix 

Two integration formulae are needed to calculate the Ps: 

/I „-2iwu 
-„^ ia = "6W (26) 

/ -2iDU 2 1 
^ (u)du= __(____,_) (27) 
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